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Abstract This article reviews the application of the 
non-equilibrium Green's function formalism to the sim- 
ulation of novel photovoltaic devices utilizing quantum 
confinement effects in low dimensional absorber struc- 
tures. It covers well-known aspects of the fundamen- 
tal NEGF theory for a system of interacting electrons, 
photons and phonons with relevance for the simulation 
of optoelectronic devices and introduces at the same 
time new approaches to the theoretical description of 
the elementary processes of photovoltaic device opera- 
tion, such as photogeneration via coherent excitonic ab- 
sorption, phonon-mediated indirect optical transitions 
or non-radiative recombination via defect states. While 
the description of the theoretical framework is kept as 
general as possible, two specific prototypical quantum 
photovoltaic devices, a single quantum well photodiode 
and a silicon-oxide based superlattice absorber, are used 
to illustrated the kind of unique insight that numerical 
simulations based on the theory are able to provide. 
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1 Introduction 

The demand for higher photovoltaic energy conversion 
efficiencies has recently led to the emergence of a whole 
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new generation of solar cell concepts based on multiple 
junctions pQ, intermediate bands [2], multiple exciton 
generation [3j or hot carrier effects [4] . Compared to the 
standard bulk technology, the corresponding devices are 
of highly increased structural and conceptual complex- 
ity, comparable to the advanced optoelectronic devices 
in the field of light emission, such as solid state lasers 
and light emitting diodes, and like these counterparts, 
they strongly rely on the design degrees of freedom 
offered by the utilization of semiconductor nanostruc- 
tures. This instance is hardly surprising if one considers 
that photovoltaics ("photons in - electrons out") can 
be regarded as the inverse regime of solid state lighting 
("electrons in - photons out"). However, even though 
the structures be the same, in switching to the photo- 
voltaic operation, the focus is shifted from the optical to 
the electronic output characteristics, with the objective 
to achieve maximum carrier extraction rate at largest 
possible separation of contact Fermi levels. 

While there is still a relation between photovoltaic 
and light-emitting quantum efficiency through the op- 
toelectronic reciprocity relation j5] , real devices are usu- 
ally non-radiatively limited, and efficient charge trans- 
port from the location of photogeneration to and fast 
extraction at the contacts becomes thus of pivotal im- 
portance. Indeed, the semiconductor nanostructures not 
only have to provide tunable and localization enhanced 
optical transitions (photogeneration) , but their absorb- 
ing states need either be coupled to extended bulk-like 
states or among themselves in order to enable electronic 
current flow (charge separation). The coupling can be 
provided via hybridization in situations where carriers 
can tunnel between localized and extended states, as in 
the case of (multi) quantum well diodes, or via inelastic 
scattering processes such as absorption and emission of 
phonons or low energy photons, respectively. The for- 
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mer mechanism is important both in devices relying on 
thermionic emission, such as multi-quantum-well solar 
cells [BJ, as well as in superlattice absorber devices in 
the sequential tunneling transport regime [7J. The lat- 
ter mechanism is utilized in the intermediate band solar 
cell, where the carriers from the isolated intermediate 
absorber band are excited into the contacted continuum 
via an additional photon absorption process. 

Apart from providing coupling mechanisms, non- 
radiative inelastic scattering mechanisms are also re- 
sponsible for any non-radiative output power losses in 
the transport process, either via recombination of charge 
carriers or via energy dissipation due to charge carrier 
thermalization and relaxation. A substantial reduction 
of this energy dissipation forms the basis of the hot 
carrier solar cell concept, which on the other hand re- 
quires the use of energy selective contacts to extract the 
hot carriers. Again, the application of semiconductor 
nanostructures is being considered, e.g. resonant tun- 
neling through single quantum dots [8]. This is an ex- 
ample of a purely electronic use of nanostructures in 
solar cells, independent of the optical properties. An- 
other example would be the interband tunnel contact 
in multi-junction devices [S]. The opposite situation of 
purely optical functionality concerns the vast field of 
light management in solar cells via photonic structures, 
which lies beyond the scope of this review. 

As a result of the large number of implications that 
are linked to the use of nanostructures in photovoltaic 
devices, it is very hard to say whether a certain struc- 
ture is going to be beneficial or detrimental to the de- 
vice characteristics, even under the idealized conditions 
of the radiative limit, assuming defect-free interfaces 
etc. In many devices, there is for instance the intrinsic 
conflict mentioned above between maximum localiza- 
tion needed for ideal absorption and maximum der- 
ealization required for ideal transport. Indeed, at this 
stage, only a few of the high-efficiency concepts have 
been realized successfully, such as the III-V semicon- 
ductor material based multi-junction devices and strain- 
balanced multi-quantum-well solar cells, the latter ex- 
ample demonstrating in an impressive way the feasibil- 
ity of highly efficient nanostructure based devices with 
a large number of interfaces. In order to be able to asses 
the potentials and capabilities of the various novel pho- 
tovoltaic concepts, a realistic theoretical estimate of the 
specific device characteristics is thus highly desirable. 

The complexity of structure and physical mecha- 
nisms as well as the prominent role of dimensional and 
quantum effects characterizing the operation of these 
novel solar cell devices preclude the use of standard 
macroscopic bulk semiconductor transport theory con- 
ventionally used in photovoltaics, which is nothing else 



than the charge continuity equation for electrons and 
holes with a source term corresponding to the net in- 
terband generation rate and a drift-diffusion current 
with a longitudinal electric field obeying Poisson's equa- 
tion. On the other hand, rate equation type quantum 
optical approaches used for solid state lasers or light- 
emitting diodes often lack a satisfactory description of 
dissipative charge transport and extraction/injection at 
contacts, while suitable quantum transport formalisms 
were devised for unipolar devices and do thus normally 
not include the optical coupling. Hence, for a com- 
prehensive theoretical description of photovoltaic de- 
vices based on nanostructured absorbers and/or con- 
ductors, a theory is required that treats on equal footing 
both quantum optics and dissipative quantum trans- 
port. These requirements can only be satisfied on the 
level of quantum kinetics, such as provided by the Wigner- 
function, density-matrix and non-equilibrium Green's 
function formalisms. 

Among the quantum-kinetic theories, the non-equilibrium 
Green's function (NEGF) method is the most versa- 
tile and powerful tool to study non-equilibrium prop- 
erties of nanostructures, since it is based on a quan- 
tum field theoretical approach to non-equilibrium sta- 
tistical mechanics [MHimUlHl] • In the NEGF ap- 
proach, the Green's functions (GFs) for the involved 
(quasi)particles (electrons, holes, phonons, photons, ex- 
citons, plasmons, etc) are the model functions provid- 
ing the physical quantities that characterize the sys- 
tem. They correspond to the response of the system 
to external perturbations, the latter entering the equa- 
tions of motion for the Green's functions, the Dyson's 
equations, in the form of self-energies. The self-energies 
due to (weak) interactions such as electron-phonon or 
electron-photon coupling are calculated using standard 
diagrammatic or operator expansion techniques used 
in many-body perturbation theory [15] . The coupling 
to the environment represented by contacts, i.e. the ap- 
plication of open boundary conditions, is provided by a 
special type of boundary self-energy. 

Due to the generality of the method steming form its 
sound foundations, the NEGF formalism has found ap- 
plication in the description of large number of different 
micro- and mesoscopic systems under non-equilibrium 
conditions, of which only a few can be mentioned here. 
Apart from the application to actual non-equilibrium 
quantum transport phenomena comprising ballistic trans- 
port and resonant tunneling in semiconductor multi- 
layers and nanostructures of different dimensionality 
(quantum wells [TBJ, wires [TTIfTS] and dots [IH]), metal- 
lic and molecular conduction [201I5I1I321I1%II241I25||2B] . phonon 
mediated inelastic and thermal transport [2"71l2lfll2"51l501 
151) . Coulomb-blockade [52H55] and Kondo-effect [Ml 
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[35] . it is also used to describe strongly non-equilibrium 
and interacting regimes in semiconductor quantum op- 
tics requiring a quantum kinetic approach [36 , 7711551 
l31fll4T)ll4"T] , with phenomena such as non-equilibrium ab- 
sorption, intcrband polarization, spontaneous emission 
and laser gain. The concept was first adapted to the 
simulation of transport in open nanoscale devices on 
the example of tunneling in metal-insulator-metal junc- 
tions 142], and has in the following been applied to 
investigation and modelling of MOSFET[43j|12|45j^ 
g7], CNT-FET [481149], resonant tunneling diodes [50] 
I5T 1 I?7 1 I5 ^ IT6 1 I53 1 and interband tunneling diodes [5411551 
[561157] . interband quantum well lasers [35] and intra- 
band quantum cascade lasers [55159], as well as infrared 
photodetectors [60], CNT-photodiodes [6411521163] and 
quantum well LEDs [54]. The theory was formulated 
both for continuum effective-mass [3"T1[45] an atomistic 
multiband tight-binding [T5] models of the electronic 
structure, and in the case of molecular conduction it is 
combined with ab-initio methods such as density func- 
tional theory [2211631124] . 

This review now provides an overview of the gen- 
eral framework for the theoretical description of quan- 
tum photovoltaic (QPV) devices based on the NEGF 
formalism, as first introduced in [661167] for the investi- 
gation of quantum well solar cells (QWSC). After the 
specification of the type of system under consideration, 
the non-equilibrium quantum statistical treatment of 
its different physical degrees of freedom (charge carri- 
ers, photons, phonons) is discussed, resulting in a gen- 
eral microscopic theory of photovoltaic devices. The 
unique kind of spectral information the theory is able to 
provide is then illustrated on the example of two proto- 
typical QPV devices, a single quantum well photodiode 
and a silicon-oxide based superlattice absorber device. 



2 Theoretical framework 

As in any device simulation, the ultimate goal of the 
present approach it to predict the response of a device 
to a given variation in the external conditions (bias, 
temperature, illumination,...) and for given material 
properties (electronic, phonemic, photonic structure), 
via the solution of (coupled) dynamical equations for 
the degrees of freedom of the device affected by the 
perturbation. In a formalized picture, the dynamical 
equations describe the evolution of the states of the sys- 
tem and their occupation for given initial solution and 
boundary conditions. Before turning to the dynamical 
equations, let us thus briefly discuss the characteristic 
set-up of the system to be described. 



intrinsic nano- 
structure absorber 




Fig. 1 Basic structure and functional elements of a generic 
quantum photovoltaic device. 



2.1 Definition of a QPV device 

Fig. [I] shows the basic building blocks of a quantum 
photovoltaic device, which are in essence the same as 
for a conventional solar cell, namely an (intrinsic) ab- 
sorber with carrier species selective contacts, usually 
obtained via selective doping resulting in a bipolar junc- 
tion device. The difference to the standard configura- 
tion is made by the presence of additional nanostructure 
absorbers in the intrinsic region, leading to the forma- 
tion of quantum wells, wires or dots and thus exhibiting 
states with increased degree of localization. The band 
diagram along the direction of transport which corre- 
sponds to such a device structure is shown schemat- 
ically in Fig. [2](a). The bulk leads are absent in the 
case where the nanostructure states form both absorber 
and conductor, i.e. in the case where the nanostruc- 
ture states are extended in transport direction, as in 
nanowires or superlattices. At the contacts, photogen- 
erated carriers are extracted and thermalized carriers 
are injected according to the splitting of the chemical 
potentials fj, ntP of electron and hole reservoirs by the 
terminal voltage Vuas- Apart from this functionality, 
which is common to all electronic devices, there is an 
aspect of contacts in photovoltaics that is even more es- 
sential, which is the mentioned carrier selectivity. It can 
be shown that even in absence of built in fields, the exis- 
tence of carrier selective contacts enables photovoltaic 
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Fig. 2 (a) Energy band diagram corresponding to a device 
similar to the structure shown in Fig. ^ The bulk leads are 
absent in the case where the nanostructure states form both 
absorber and conductor. The device (absorber/emitter) states 
can be bound, quasi-bound or form a quasi-continuum. At 
the contacts, photogenerated carriers are extracted and ther- 
malized carriers are injected according to the chemical po- 
tentials ii n ,p for electrons and holes, respectively. In the de- 
vice regions, the physical processes relevant for the photo- 
voltaic operation are: (b) radiative transitions, i.e. (1) inter- 
band photogeneration and radiative recombination as well as 
(2) radiative intraband transitions, which may also lead to (3) 
photon recycling; (c) coherent and dissipative quantum trans- 
port involving non-radiative intraband transitions, such as (1) 
phonon mediated carrier capture and escape, (2) intraband 
relaxation, (3) scattering assisted or (4) direct tunneling be- 
tween absorber states, and (5) phonon assisted carrier escape 
as combination of the two; (d) non-radiative recombination, 
via (1) the Auger mechanism or (2) deep trap states. 



device operation [68] . For some photovoltaic applica- 
tions like the hot carrier solar cell, contacts need in ad- 
dition be energy selective. Even though most solar cells 
are rather macroscopic than mesoscopic devices and the 
transport properties are not primarily determined by 
the contact resistance, the contacts are usually criti- 
cal components in the solar cell design. This is due to 
electronic matching problems, leading to band offsets, 
Schottky barriers, strong local fields and defect forma- 
tion. Indeed, defect recombination at contacts repre- 
sents a major loss mechanism, and large efforts are 



dedicated to the passivation of these surface defects. 
While the contact states do normally form a continuum, 
the absorber states in the device region can range from 
completely localized bound states to quasibound states 
or even extended quasicontinuum states, depending on 
the degree of confinement imposed by the nanostruc- 
ture, and they enable a large number of additional ra- 
diative and non-radiative transitions which are relevant 
for the photovoltaic device operation, as shown in Figs. 
[2p3)-(d). Here, the transitions which do involve only one 
carrier species and conserve particle number are termed 
"intraband" processes, even though they take place be- 
tween discrete states rather within a proper band or 
even between subbands. Analogously, transitions across 
the fundamental gap ar termed "interband" transitions. 
In the case of absorber states with reduced extension 
in transport direction, non-local processes such as di- 
rect or scattering assisted tunneling connecting spa- 
tially separated states or localized to extended states 
become pivotal for the device performance, and must 
thus be considered in any appropriate simulation of 
such a device. 



2.2 NEGF description of a QPV system 

Due to the large number of applications of the NEGF 
formalism, its foundations have been discussed in many 
reviews and books, e.g. [69,70,71,40], wherefore we will 
limit the discussion here to the aspects that are relevant 
for photovoltaic applications. Many of these aspects are 
also present in resonant-tunneling devices on the quan- 
tum transport side and in light-emitting devices on the 
quantum optics side, and the discussion thus follows 
some of the lines of the corresponding literature, such 
as presented in [16] and [38]. We start the discussion 
with the quantum-statistical mechanics picture of the 
semiconductor device. 

2.2.1 Field operators and Green's functions 

A comprehensive, unified description of the quantum 
photovoltaic device is found on the quantum-statistical 
mechanics level, where the optical, electronic and vibra- 
tional degrees of freedom of the system are described by 
the corresponding quantum fields for photons, charge 
carriers and phonons. Within the NEGF theory of quan- 
tum optics and transport in excited semiconductor nanos- 
tructures, physical quantities are expressed in terms of 
quantum statistical ensemble averages of single parti- 
cle operators for the interacting quasiparticles intro- 
duced above, namely the fermion field operator for 
the charge carriers, the quantized vector potential A 
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for the transverse photons and the ionic displacement 
field U for the phonons. The corresponding GFs are 



G(l,2) 
2>L(1,2) 



h 

1 i 
Mo ft 



{A,{l)A v {2)) c 
- {A,{l)){A v {2)) 



(electrons) (1) 



(photons) (2) 
(phonons) (3) 



where (...)c denotes the contour ordered operator av- 
erage peculiar to non-equilibrium quantum statistical 
mechanics [T2|[T3] for arguments 1 = (ri,^) with tj 
on the Keldysh contour [13]. In the case of the lat- 
tice displacement field, the continuous spatial coordi- 
nate r is replaced by the discrete ion position vectors 
R = L + k = Lk, where L is the equilibrium lattice 
vector and k, the basis vector within the unit cell [ID] . 
The associated four-vector is 1 = (LiKi,]^). 

Since it is a declared goal of the present approach to 
treat systems where any of the degrees of freedom may 
be quantized to some extent via dimensional confine- 
ment, we will keep the treatment as general as possible 
and thus use the real-space, complex-time representa- 
tion of the operators and their averages in the review 
of the formalism, even though in many cases the choice 
of a reciprocal or modal space reduces formal and com- 
putational complexity significantly. 

2.2.2 Dynamical equations 

The GFs are determined as the solutions of correspond- 
ing Dyson's equations [7"2"ll3"gHU] . 

f d3 [Go 1 (1, 3) - E(l, 3)] G(3, 2) = 6(1 - 2), 

f rf3 [CPSr'ai) - nP(l,3)] ^(3,2) = 5(1-2), 
Jc 

(4) 

where the integration is J c dl = J c dt\ J d 3 ri. The GFs 
Go, T>q and T>q are the propagators for non-interacting 
electrons, photons and phonons, respectively «-» de- 
notes transverse and boldface tensorial quantities. The 
electronic self-energy £ encodes the renormalization of 
the charge carrier GFs due to the interactions with pho- 
tons and phonons and other carriers, and enables thus 
the description of charge carrier generation, recombina- 
tion and relaxation processes. It is responsible for the 
appearance of excitonic effects in the carrier spectrum 



^(3,2) = V(l-2) 



and leads to band-gap renormalization under high ex- 
citation. Charge injection and extraction at contacts is 
considered via an additional boundary self-energy term 
reflecting the openness of the system. The photon and 
phonon self-energy tensors and JJ V describe the 
renormalization of the optical and vibrational excita- 
tion modes due to the interaction with the electronic 
system, i.e. absorption and emission of photons and 
phonons, leading to phenomena such as photon recy- 
cling or hot carrier effects and including excitonic sig- 
natures in the bosonic spectra. 

The self-energies can be derived either via perturba- 
tive methods using a diagrammatic approach or a Wick 
factorization, or using variational derivatives. Their de- 
pendence on the GFs can be expressed using Hedin's 
approach [73 (henceforth, contour integration over in- 
ternal variables is assumed) 

E(l, 2) = - ihQG(l, 3)r e (3, 2, 4) W ef f (4, 1) 

- ihj(l, l')G(l, 3)F? (3, 2, 4) V(4, 1011=1', 

(5) 

2) =ihj(l, l')G(l, 3)1^ (3, 4, 2)G(4, F)| 1=1 , , 

(6) 

n p (l2) =^a,l)G(l,3)r p (3,4,2)G(4,l), (7) 

where Q, and T are the bare interaction vertices. In 
the expression for the phonon self-energy, the RHS is 
integrated over all continuous spatial coordinates. The 
effective carrier-carrier interaction potential W e / / is the 
sum of an electronic plasmon screening and a phonon- 
mediated contribution, 

W ef f ( 1 , 2) = We (L 2) + W p (1,2), (8) 
Wed, 2) = W e(0) (I, 2) + W e(0)) (I, 4>(4, 3)W e (3, 2), 

(9) 

W p (l,2) = W e (l,4)E> n (3,4)We(4,2). (10) 

Here, 7r is the longitudinal electromagnetic polarization 
function corresponding to the self-energy of longitudi- 
nal photons and obeying the equation 



tt(1,2) = ih QG(l,3)r e (3,4,2)G(4, 1) 



(11) 



and T> n is the density-density correlation function of 
the nuclei, which is related to the phonon GF in real 
space via [I0J 

2?"M,iV) = E E V Q p K (r-L-«) 

x V p ap (Ln, L'k', t, t')V/3p K >(r' - L' - «')> (12) 
where p K is the local charge density due to the ion cores. 
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The equations for the vertex corrections used in the 
iterative determination of the self-energies are 

r e (I,2,3) = - Q5(1,2)<5(1,3) + ^k^G(4,6) 

xr«(6,7,3)G(7,5), (13) 
r 7 (I,2,3) = - VoJ(hl')d(Lmh3)\v=i 
, ^(1,2) 



<5G(4,5) 



G(4,6)r^(6,7,3)G(7,5), (14) 



r*(l,2,3) = - ^(1)5(1,2)5(1,3) + g|l||G(4, 6) 

x r p (6,7,3)G(7,5). (15) 

Inserting these expressions in the equations for the po- 
larization functions provides self-consistent approxima- 
tions for the latter in the form of the Bethe-Salpeter 
equation for the scalar four-point polarization function 

V{1, 1', 2, 2') =7yi, I', 2, 2') + V Q (L 3, 4, 1') 

x ^(3,4,5,6)^(5,6,2,2'), (16) 
511(3, 4) 



^^-*G(5,6r 
through the relations 



17) 



vr(l, 2) =ihQ 2 V(l, I, 2, 2), (18) 
^ 7 (I,2) =^a.3)J r (4 ! 4 / )P(l,3,2,4)| 4 , =4 , (19) 



^ P (L I) =ihT(X, 1)^(2, 2)-p(l, 1, 2, 2). 



(20) 



For practical evaluation, the above equations are 
rewritten in terms of the real time components of the 
GFs and self-energies according to the order of the time 
arguments on the contour [T51IT4"] . i.e. the lesser (<), 
greater (>), retarded (R) and advanced (A) compo- 
nents. In photovoltaics, in contrast to standard appli- 
cations of NEGF in ultra-fast laser physics, but similar 
to many other quantum transport situations, one is pri- 
marily interested in the steady-state behaviouiQ and it 
is thus usually sufficient to consider the dependence on 
the time difference r = t\ — t 2 , which permits a trans- 
lation of the investigation into the energy domain via 
the Fourier transform 

/CO 
dTei ET O a (r 1 ,r 2 ; T ), (21) 
-OO 

O e {(G, D 7,p ), (17, IP*)}, a € {R,A, <, >}, 

1 There are photovoltaic devices where the investigation of 
time-dependent processes is crucial, e.g. the ultra-fast injec- 
tion of charge from the dye into the semiconductor in dye- 
sensitized solar cells. 



2.3 Hamiltonian and self-energies 
2.3.1 System partitioning 

In the formulation of the system Hamiltonian, a com- 
mon approach is to partition the system into easily ac- 
cessible term, with most of the complexity in the re- 
mainder. In the present case, there are three aspects 
of complexity that may be tackled by this way, namely 

1. the openness of the system, i.e. the requirement to 
include the contacts in the description of the system, 

2. the presence of interactions in the device, and 3. the 
different degrees of freedoms to be considered, i.e. elec- 
tronic, optical and vibrational. The first aspect leads to 
a system Hamiltonian partitioning of the form 



H =H L 



H H + H 



DR 



(22) 



consisting of the separate terms for isolated device, con- 
tact reservoir and the coupling between the two. In the 
second case, the system Hamiltonian is written as the 
sum of non-interacting and interaction terms, 



H =Hq + Hi 



(23) 



Finally, we have the different components attributed to 
electrons, photons and phonons, 



H =H e + Hy + Hp, 



(24) 



where the terms include the mutual interactions as well 
as the interactions with further system components such 
as impurities etc. 

In order to make the system treatable, interactions 
in the contacts are not considered explicitly, and their 
effect on the contact-device coupling is neglected, which 
amounts to the assumption Hf — Hf R = 0. Further- 
more, the bosonic subsystems will be regarded as non- 
interacting, i.e. i/i i77 = H i pp = 0. This leaves the total 
system Hamiltonian 



H 
H a 



--H + Hi 



■Hi 



H, 



DR 



HP =H° 



H. 



D 



H, 



D 



(25) 
(26) 

(27) 



To quantify the action of the Hamiltonian on the 
many-body system, it is represented in second quanti- 
zation using the field operators introduced in Sec |2.2.i] 
e.g. for the electronic system 



H{t) = J d 3 r&(r,t)H e <P(r,t). 



(28) 



In the case of the carrier-carrier interaction, this pro- 
duces a two-particle Hamiltonian, while all the other 
terms are on the single-particle level. 
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The Hamiltonian can now be used in the derivation 
of the equations for the self-energies using the many- 
body perturbation theoretical expansion of 



G(r,t;r',t') = -~ (e~i So <^'M£(r, i)<pt(r'> *')) 



and 



(29) 



-R/c rf ^'( s )jr«( r ^)j-«t( r ' jt ') 
(Mo 1 ,!), F'" = (AW), (30) 



In photovoltaic devices, the electronic system is open 
in at least one dimension (r'), for which thus scattering 
states have to be used rather than eigenstates of the 
closed system. However, the remaining confined dimen- 
sions (f) may still be described by the corresponding 
eigensolutions, allowing for an eigenmode expansion of 
Hamiltonian, GFs and self-energies in these dimensions, 
which results in a general field operator representation 



r = (f , r'). 



(35) 



where the perturbation is given by %' = W—Wq. By re- 
stricting the perturbation to specific parts of the Hamil- 
tonian, the corresponding self-energy terms can be iden- 
tified in the resulting Dyson equation. In the following, 
the Hamiltonian terms and corresponding self-energies 
for the contacts and the various interactions shall be 
discussed. 

2.3.2 Non-interacting isolated subsystems 

Electrons The electronic system without interactions]^] 
between valence electrons and coupling to the bosonic 
degrees of freedom is described by 



photons While absorption in most nanostructured solar 
cells may be reasonably described by the linear response 
to a coherent excitation, spontaneous emission as the 
most fundamental loss process relies on the coupling to 
the incoherent field due to vacuum fluctuations. In gen- 
eral, the total electromagnetic field may thus be written 
as the superposition A = A co ^+Ai„ c of a coherent con- 
tribution A-coh = (A) corresponding to coherent light 
sources and a contribution Ai„ c with (Aj„ c ) = from 
incoherent light sources and from spontaneous emis- 
sion. While the coherent vector potential is related to 
the time-dependent part of the classical electric field £ 
via the standard relation 



2m 



(31) 



where mo is the free electron mass. In the above ex- 
pression, the first term provides the kinetic energy and 
U contains the potential for the interaction of valence 
electrons with the ion cores. It is further customary to 
include in U also the Hartree term of the Coulomb inter- 
action corresponding to the solution of Poisson's equa- 
tion that considers carrier-carrier interactions (H ee ) on 
a mean-field level. The GFs corresponding to this Hamil- 
tonian is defined via 







G (l,l') = 5(1,1'), 
with the inverse function in Q given by 



G - 1 (l ! l') = 



5(1,1')- 



(32) 



(33) 



The Hamiltonian operator in ( 31 1 is the one entering 



the steady-state Schrodinger equation corresponding to 
the single particle eigenvalue problem of the isolated 
system, 



(34) 



2 In contrast to the case of photons and phonons, there 
are no completely free electrons in the system, since there is 
always the lattice potential due to the ion cores. 



~A co/t (r,i) 



(36) 



and is thus obtained from the solutions of Maxwell's 
equations according to classical electrodynamics, the 
incoherent part needs to be treated on the quantum 
mechanical level via the GF introduced above and will 
be discussed in what follows. 

In terms of the vector potential operator introduced 
at the beginning of this section, the Hamiltonian of 
transverse electromagnetic radiation reads |40j 



(37)" 



The free contour-ordered photon GF in real space obeys 
the equation 



1 d 2 



#o(U') = V(i,i'). 



(38) 



which, in analogy to (33), defines the inverse free pho- 



ton GF entering the photon Dyson equation in Q via 
1 d 2 



^(M') 



c 2 d% 



V(i,i') 



(39) 
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Since it is a solution to the wave equation following 
from Maxwell's equations, the quantized field can be 
represented as the free-field mode expansion 



For normal incidence along r' (f _L f) and thermal 
radiation from source at temperature T, the occupation 
can be written [7U 



A(r, t) = £ [M\ q)aA,q(*) + AS(A, 



-q)4,- q (*) 



A.q 



^A,q = n BE (fiCo\q\,T)e(q , )8(q'ta,na- |q| 



(49) 



x e 



qr 



A (A,q) = /o „ fc £ Aq, 



(40) 
(41) 



where e Aq is the polarization of the photon with wave 
vector q and energy fiw q added to or removed from 
photon mode (A, q) by the bosonic creation and anni- 
hilation operators 



where tlbe is the Bose-Einstein distribution function 
and a is half of the aperture angle under which the 
source is seen by the solar cell. 

In the general case of an inhomogeneous medium, 
one may opt to include the modification of the opti- 
cal mode due to spatially varying background dielectric 



(42) 



and V is the volume. In this representation, the non- 
interacting Hamiltonian of transverse radiation acquires 
the simple form 



constant £&(r) in (39) for the "free" propagator by re- 
placing Co by c(r) = co/n Q (r), where n = ^Jei- In the 
waveguide situation, the resulting problem gives rise to 
confined optical modes, which, in analogy to the elec- 
tronic case, can be used to expand the photon GFs [75] , 
As discussed above, solar cells are optically open in at 
least one dimension (r r ). The general expression for the 
photon field operator is thus 



A,q 



1 



(43) 



A(r,i) = ^^(f)A n (r',t), r=(f,r'), 



(50) 



With the representation Q, the EM GF in ^ can be 
related to the bare scalar photon GF corresponding to 
the following expectation value of bosonic operators, 



7,(0) 



V 



via 



(q;t,t') = - % - (a A>q (t)a A>q (t') + a A _ q (t)a A ,_ q (t') 

(44) 



2e V"wA, q 



e fc Aq fc Aq 



Noting that the photon occupation number operator for 
free field mode (A, q) is given by n Aq = a Aq aA, q , the 
components of the scalar photon GF can be expressed 
in the energy domain as 



where V?n( r ) are the eigenmodes in the confined dimen- 
sions. 

phonons In analogy to the photon field, the lattice dis- 
placement field for a bulk semiconductor with full trans- 
clational symmetry can be expanded into vibrational 
eigenmodes in reciprocal space, 

U a {LK, t) = ££W 0) a«(A Q) [h,Q(t) + b\ t _ Q {t)] 
Q A 

x e ,Q ' (L+K) , (51) 

h 

y/2NM K hn A a' 

(a = x,y,z), 



( 45 ) U 0<aK (A,Q) = 



i(Q), 



(52) 



2^>(q; E)=- 2-ki iq 5(E T RwA.q) 



V 



i,R/A 



{Nl^ + l)5{E±hw x . CL ) 
1 



(46) 



where N is the number of atoms, M K the mass of the 
basis atom at k, and ^a.Qj^q arc the bosonic cre- 
ation and annihilation operators for a phonon mode 
(A, Q) with polarization e QKy i(Q) and energy HHa.q 
determining the time-evolution in complete analogy to 



( 42 ) . The phonon modes of a given structure follow from 



with JV£ q 



E — huj\ q d= Z77 E + hu>x qilf) 

(47) 

the mode occupation, which in the 



'A,q/ 



free-field case is related to the incident photon flux via 



the eigenvalue equation 

£X QK:/ 3 K '(Q)ea K /i(Q) =(^i,Q) 2 e QK A(Q), (53) 

in terms of the dynamical matrix Ai, which in turn is 
determined by the interatomic force constants $ via 



K« = US/Co. 



(48) 



The modal photon flux in turn is given by the modal 
intensity of the EM field through A = 7 A q /(/twx q ). 



M 



(Q) 



\ZMkM; 



_^ e -iQ-(L+«-«')^(L K)0 Av). 

(54) 
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The interatomic force constants can be obtained via 
phenomenological models or ab-initio calculations |76[ 
177) . Again, it is possible to express the Hamiltonian and 
GF of the displacement field given in ^ in terms of the 
vibrational eigenmodes, the former acquiring the form 



(55) 



and in the case of the latter via the bare scalar bosonic 
propagator 



2^(Q;t,t') 



&A,Q(t)&i Q (t')+&; 



- Q (t)^,- Q (t'; 

(56) 



through the relation 



v p . ,(i,i')=y 



2e N 2 M K M K ,f2 A , Q 



xe a ^(Q)e^ M (Q)^(Q;t 1 »ti')- (57) 

In the case of coupling to equilibrium phonons, the 
scalar propagator can be expressed in a form similar 



to Eq. (47), with the occupation q given by the 
Bose-Einstein distribution at the temperature of the 
heat bath. 

In the case of broken translational invariance caused 
by the use of nanostructures, the size of the eigenvalue 



problem in ( 53 1 will increase by the dimension of the 



unit cell in the direction of aperiodicity, while the di- 
mension of the reciprocal space is correspondingly re- 
duced. However, the system is not elastically open in 
the sense that a scattering problem needs to be solved, 
i.e., the vibrational state can be described in terms of 
vibrational eigenmodes in all dimensions. 

2.3.3 Contacts 

In theory, contacts represent the open boundary con- 
ditions that need to be imposed on an open system in 
contact with the environment in order to enable charge 
exchange between the device and the reservoir formed 
by the contact. In macroscopic theories, boundary con- 
ditions are applied to the carrier density: surface charge 
is either zero in the case of ohmic contacts (local charge 
neutrality) or fixed by the Schottky barrier height in the 
case of metal-semiconductor contacts. In microscopic 
transport theories, the description of contacts to open 
systems is itself a whole field of ongoing research. At 
the most basic level, the formalism is required to pro- 
vide information about the density of states in the con- 
tact, the occupation of these states, and the coupling 
of the contact states to the states in the device. In the 



NEGF formalism, all this information is contained in 
the boundary or contact self-energy. 

In the electronic case, the general treatment (see 
e.g. [7T]) starts from the partitioning of the overall GFs 
into contact and device contributions, according to the 



overall Hamiltonian in ( 22 ) (in a discrete basis) 



G = 



(E + irj)l-~H. D -t 

(E + iri)l -H fl . 



-rt 



DR. 

r<R n R 

^RD ^R 



(58) 
(59) 



where r = Hdr- Eq. (59) is identical to the Dyson 
c equation produced by a first order perturbation in the 
contact-device Hamiltonian, i.e. using T-L' = Hdr in 



(29). If written out in components, it can be used to 



obtain the retarded device GF G R in terms of the re- 
tarded GF 



g R = [(E + itfl-Hn] 



(60) 



of the isolated reservoir: replacing G R R = G^rg^, the 
retarded device GF results in 



Gg = [(£ + ^ ? 7)l-H D -T+g fi T] (61) 

which defines the retarded boundary self-energy 

V RB = t^ h t. (62) 

With the assumption of vanishing scattering between 
contacts and device, the Keldysh relation for the cor- 
relation functions of device and reservoirs, in analogy 



to ( 59 ) provides the in- and outscattering contact self- 



energy terms 



^ t 
te>t\ 



(63) 



where g§= is the lesser/greater GF of the uncoupled 
reservoir. A frequently used assumption is that of con- 
tacts which are equilibrated due to scattering, which is 
usually a good approximation for solar cell devices oper- 
ating at elevated temperature and possessing extended 
contacts. In this case, the lead correlation functions can 
be expressed in terms of the retarded and advanced 
functions using the fluctuation-dissipation theorem, 



g< =iU(E)a, g> = i[l - U(E)]a, 



(64) 



where /i is the chemical potential and a = i[g R — g ] 
the spectral function of the reservoir. 

The whole problem of the coupling to the reservoirs 
is now reduced to the calculation of the contact GF 
g fl , which in principle is of infinite dimension, but only 
needs to be known in the close vicinity of the device 
boundary, owing to the reduced dimensionality of the 
coupling matrix r, and can therefore be calculated by 
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surface GF methods using decimation techniques [78] 
[79] , conformal maps [80] or complex band methods [81] 
15211551155] . 

Similar surface GF approaches can also be used for 
photons [75,84 or phonons [551I5B] in systems that are 
open in the optical or vibrational sense, respectively. 

2.3.4 Interactions 

The effects of interactions between the different degrees 
of freedom can be considered perturbatively via respec- 
tive self-energy terms in the same way as the pertur- 
bation due to contacts by replacing the perturbation 



Hamiltonian in (29) and (30) with the corresponding 



interaction term. Alternatively, the self-energies may be 
derived using variational-derivative techniques [121140) . 

electron-photon The single particle potential for linear 
coupling of the electromagnetic field to the electronic 
system is obtained in Coulomb-gauge as 



H e Jl) 



A(l)-p(l) 



(65) 



with A the vector potential operator of the electromag- 
netic field and p(l) = —ihW ri the momentum operator. 
To first order in the vector potential, the interaction re- 
sults in the singular self-energy term 



ihe 
too 



A(l)).p(l) 



ihe 
mo 



A coh (l)-p(l) 



(66) 



describing the coupling to classical light, which corre- 
sponds to a mean-field treatment of the interaction. 
Since this self-energy is proportional to the coherent 
polarization, it leads to off-diagonal entries if a band 
basis is used [55] , which can be eliminated using a band- 
decoupling scheme [72], resulting in an effective inter- 
band self-energy that is second-order in the field [?], 
e.g. for a two-band model (b — c, v) 



2 e j,cc(l, l0 — 2 & e ~ CV {V)G VV {1, lO-^Wucd')- 



(67) 



To second order in the vector potential, the self-energy 
from the Fock term of the self-consistent Born approx- 
imation is 



ih[i§e 2 ■sr-^ 



[p a {l) - p a {l')}G{\,2) 



x pp(2)^ a p% l0|i'=i, 



where the photon GF is defined in ([2]). Close inspection 
of the two expressions shows that the self-energy com- 



In analogy to the formal treatment of the electronic 
system, a perturbative or variational approach can be 
used to derive the photon self-energy component due to 
electron-photon interaction, according to the expansion 



in (30). To lowest order, corresponding to the random- 



phase-approximation (RPA) for the polarization func- 
tion, the transverse photon self-energy reads 



^(1,2) = 



ih[i§e 2 
4toq 



[p a (l)-p a (l')]G(l,2) 



x^(2)G(2,l')|i 



(69) 



Electron- Phonon The effects of the vibrational modes 
on the electronic degrees of freedom of the system are 
described in terms of the coupling of the force field of 
the electron-ion potential V e i to the quantized field U 
of the ionic displacement |87j , 



H ep (r, i) = £tf(L + k, t) • W ei [r - (L + «)], (70) 
with the displacement field given by the normal-mode 



expansion (51). The potential felt by electrons in het- 



erostructure states due to coupling to bulk phonons can 
thus be written as 

H ep (r,t) = J2 %^ Q ' r {6A,Q« + ^,_ Q (*)}, (71) 

AQ 



where r is the electron coordinate, and Ua.q are re- 
lated to the Fourier coefficients of the electron-ion po- 
tential V e i [H7]. The electron-phonon self-energy can 



be obtained from the interaction potential ( 70 ) within 



many-body perturbation theory in complete analogy to 
the electron-photon interaction. The self-energy corre- 
sponding to the Fock term of the SCBA appears in the 
form 

2*0, 2.)=J2J2J2^^- Lk ) ^ ^ - L v ) 

x^(L K ,t i; LV,t 2 )G(l,2), (72) 

with T = W e i. For inhomogeneous systems, one has 
to consider also the Hartree term [55], which reads 

2)=[EEE / ( r i - Lk ) - LV 

L,k. L' ,tz' a./3 

x D^(L K ,t 1 ;LV,t 3 )G(3,3)]i(i,2) ! (73) 

In the case where the phonon GFs are not computed 
self-consistently, it is often more convenient to use in 



ponent in ( 68 ) due to the subtraction of the coherent the derivation of the self-energy the form ( 71 1 for the in 



part is formally identical to (67). 



teraction potential. Finally, the RPA-expression of the 
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phonon self-energy required for the renormalization of 
the phonon propagator reads 



m v (1,2) =ih / d 3 n / d 3 r 2 ^ (n - Liki 



(74) 



Carrier- carrier In general, there are two types of elec- 
tronic interactions to be considered in solar cells: electron- 
electron and electron-hole interactions. While the for- 
mer affect primarily the intraband-dynamics and are 
thus of relevance in the discussion of thermalization 
and relaxation issues in transport, the latter can have 
a strong impact on the optical properties via the ex- 
citonic enhancement of interband transitions. Since we 
do not consider the case of excitonic solar cells where 
exciton diffusion is a key process, like in organic photo- 
voltaics, this absorption enhancement will be the only 
excitonic effect to be discussed. 

The most important effect of electronic interactions 
in junction-based inorganic solar cells is the built-in po- 
tential due to bipolar doping and the modification of 
the electrostatic potential U under illumination. Both 
cases can be treated macroscopically via the solution of 
Poisson's equation 



e V r [e(r)V P tf(r)] = n(r) - N d {v), 



(75) 



for a given profile Nj, of fully ionized dopands and 
(static) longitudinal dielectric function e and with the 
charge carrier density n (including holes) from the GFs. 
In the microscopic picture, the result of this approach 
corresponds to the Hartree-term, i.e. a mean-field level 
treatment, of carrier-carrier interaction. The correspond- 
ing effective single-particle interaction potential is 



Vee(r,t)= / d 3 r'&(r',t)V ee (r>-r)£(r>,t), 



(76) 



where V ee (r' — r) is the two-body interaction potential. 
The Hartree self-energy for this potential is 



Eg(l, 2) =iM(l,2) j d2K e (n - r 2 )G(2, 2) 



^(1,2)C/(1), 



(77) 
(78) 



where U is the mean-field potential. On the same level 
of perturbation theory, the corresponding Fock term 
reads 



^f e (1, 2) =iM(tj - t 2 )Ke(ri - r 2 )G(I, 2). 



(79) 



Electron-electron interactions beyond the mean-field 
Hartree-Fock level may be considered by explicit inclu- 
sion of dynamical screening via the GW formalism. For 
practical purposes, often a statically screened poten- 
tial is assumed together with the above forms of the 



self-energy. Concerning the electron-hole interaction be- 
yond Poisson's equation, excitonic effects in semicon- 
ductor nanostructures have been discussed for steady 
state linear absorption or in the regimes of high pulse 
excitations. For the latter, sophisticated quantum-kinetic 
theories were developed [5^1l371l411[72l[9Ul[75] . For extended 
systems, the computational cost of such approaches lim- 
its their applicability to special cases of reduced com- 
plexity, such as systems that may be described with a 
few eigenstates, and they are thus not suited for the 
investigation of transport properties. For this reason, 
excitonic effects will be considered here only in terms 
of the coulomb-enhancement of the optical interband 
self-energy due to the coherent polarization of the semi- 
conductor system resulting from the coupling to a co- 



herent (classical) light field, as introduced in (67). The 
singular interband Coulomb term is (a ^ b = c, v) 



S s ee . 6 (n, r 2) t) = thVir, - r 2 )G< (n, r 2 ,t, 



(80) 



where V is the (screened) Coulomb potential, and de- 
pends thus on the coherent interband polarization through 
the interband GF, which the decoupling provides in the 
form 



G vv (1, 2) S vc (2 , 3) G^ c (3, 1/ ) 



GUht)= d2 di 
Jc Jc 

+ G< (1,2)2&(2, 3) Gt (3,10 



i 

h 



(81) 

J^d2^d3 S s vc % 3)T, C (1, 2, 3, 1'), (82) 

with the tilde denoting the GF without interband cou- 
pling. Inserting the explicit expressions for the singu- 
lar self-energies leads to the Bethe-Salpeter equation 
(BSE) for the coherent polarization, which for steady 
state reads 

G< c (r 1 ,r 1 ,;E) = G< c{0) (r 1 ,r 1 ,;E) 

+ J d 3 r 2 J d 3 r 3 V(r 2 -r 3 ) 

x% c (ri,r 2 , r 3 , r v ;E)G< c (r 2 , r 3 , E), (83) 



with 

T vc,(0) 



G±^(r u r v ;E) 



. d 3 r 2 k coh {v 2 ,E) ■ p(r 2 ) 
mofi J 

X Tvc 

(ri, r 2 , r 2 , T\>] E) (84) 



the coherent interband polarization of non-interacting 
electron-hole pairs. The corresponding self-consistent 
BSE-type equation for the singular self-energy is ob- 



tained by using (83) and (84) in Eq. (80) 



S^rur'^E) = S^r^E) + Vfa - r[) J 
x J d 3 r 3 r TO (r 1 ,r 2 ,r 3 ,r 1 ,;£;)Z , f ;) (r 2 ,r 3 ,£;), 



d 3 r 2 



(85) 
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where S s ab {Q) 



S 



2.4 Steady state device characteristics from NEGF 

Due to the very general grounds of their determina- 
tion via quantum kinetic equations, the NEGF contain 
information not only about the electronic, optical and 
vibrational states of the system, but also about the oc- 
cupation of these states and the rates of transitions be- 
tween them. Knowledge of the GFs for carriers, pho- 
tons and phonons thus immediately provides informa- 
tion about the device characteristics as they are related 
directly to physical observables: the ensemble average 
of any single-body operator O can be written as [TS] 

(d(r,t))=Tr[p6(r,t)} 

= =Flim lim Tr[(9(r, t) J"<(r, t; r', t% (86) 

where T = G, T> 1,p and the upper (lower) sign applies 
to fermions (bosons). 

In photovoltaics, the properties of interest deter- 
mined experimentally cover on the one hand a range 
of spectral quantities, such as absorption, emission and 
transmission spectra, or spectral response and quan- 
tum efficiency. These quantities can be obtained from 
the GF in the energy (frequency) domain. On the other 
hand, there are the integral quantities which follow from 
the corresponding spectral versions upon energy inte- 
gration, such as the charge carrier current density pro- 
viding the photovoltaic current-voltage (J/V) charac- 
teristics which determine the energy-conversion efficiency 
of the solar cell device. 

2.4-1 Density and current 

Following the above procedure, the electron density in 
steady-state is computed by inserting the electron den- 



sity operator p(r,t) = Srt(r,i)#(r,t) in Eq. Q and 
performing the Fourier-transform to the energy domain, 



n e (r)=-i J ^G<(t,t;E). 



(87) 



The hole density is obtained by replacing in the above 
expression G < with —G > . In the same way, the photon 
correlation function yields the photon density via 



(r) 



2tt 



Tr{2r<(r,r,£:)}, 



(88) 



where the integrand is given by Planck's law if the cor- 
relation function for equilibrium free-field photons is 
used. 



The electronic current density is determined as the 
ensemble average of the operator 



je(r,*) 



-ieh r 



2mo 
= lim 



& (r)V r !^(r) - {Vr# f (r)}#(r) (89) 
ieh 



2 mo 



[V r - V r #*(r)!f(r'), 



(90) 



and in steady state can be expressed in terms of the 
carrier GF via 

je(r)= lim JL^-VJ f ~G<(r,r';E), (91) 

This form of the current may also be obtained directly 
from the formulation of the electron continuity equa- 
tion in terms of the electronic GFs by inspection of the 
divergence term. 

The photon current j 7 in the device follows from 
the evolution of the light intensity, which is given by 
the symmetrized non-equilibrium ensemble average of 
the Poynting vector operator |91j . 



l(S(rx))| 
where 



d(huj)huj j 7 (r 1; hui), 



1 I- 



(S(r!,ti)) =— (£(ri,ii) x8(r 2 ,t 2 ) 
Mo \ 

with 



sym 1 l=2 



g (r>t) = _a^ 1 t) ) B( r ,t)=Vx A(r,<) 



(92) 



(93) 



(94) 



the electric and magnetic components of the EM field. 
Recalling the definition Q of the photon GF, the spec- 
tral photon current results in terms of the latter to 



3i,i 



(ri,fcd)= lim — V Re[v,(2){l?> (n,r 2 ;^) 

+ 2?< (r x , r 2 ; fej)} - V„ (2) {D> v {v x , r 2 ; foj) 



^(ri,r 2 ;?iw)} 



2-4-2 Conservation laws and scattering rates 



(95) 



The macroscopic balance equation for a photovoltaic 
system is the steady state continuity equation for the 
charge carrier density 



V-jc(r) =Q c (t) -TZ c (r), c = e,h, 



(96) 



where j c is the particle current density, Q c the (volume) 
generation rate and 1Z C the recombination rate of carri- 
ers species c . In the microscopic theory, the divergence 
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of the (particle) current density is given by [12|fT3] 



Vj(r 



dV 



E R (r,r';E)G < (r',r;E) 



2 r dE 
~VJ 2nh, 
■E<(r, r'; E)G A {r' , r; E) - G R (r, r'; E) 



xE < (r',r;E)-G < (r,r';E)E A (r , ,r;E) . (97) 

The RHS of the above equation provides a general ex- 
pression for the microscopic intra- and interband rates 
of (potentially non-local) scattering processes. Current 
conservation requires the RHS to vanish for energy in- 
tegration over all bands or states, which means that 
any scattering contributions originating in transitions 
in energy space must cancel upon energy integration. If 
the integration is restricted to either conduction or va- 
lence bands, the above equation corresponds to the mi- 



croscopic version of ( 96 ) and provides on the RHS the 



total local interband scattering rate. The total inter- 
band current is found by integrating the carrier specific 
divergence over the active volume, and is equivalent to 
the total global transition rate and, via the Gauss the- 
orem, to the difference of the interband currents at the 
boundaries of the interacting region. Making use of the 
cyclic property of the trace, it can be expressed in the 
form 



TZfnt — 



d 3 r I I d 3 r' 



dE 

V J ' I 27rk 

Z > (r,r';E)G < (r',r; E) 



E<(r > r';E)G>(v',r;E) 



(98) 



with units [TZ] — s" 1 . The above equation can be used 
to define the spectral scattering current J<p(E) via IZ^, = 
J dEJ ! j ) (E) for a given scattering process described by 
a self-energy E^, such that [7T] 



ME) = 



d 3 r I d 3 r' 



1 

-E>(r,r';E)G<(r',r;E) 



S<(r,r';E)G>(r',r;E) 



(99) 



= ;r-rtr 
2ttH 



E<(E)G>(E)~E>(E)G<(E) , (100) 



which describes the total scattering rate at a given en- 
ergy as difference of the inscattering rate into an empty 
state at a that energy h~ 1 tv\S^{E)G > {E)\ and the 
outscattering rate fi. _1 tr[X'f : (£')G > (£')] from an occu- 
pied state at the same energy. 

If we are interested in the interband scattering rate, 



we can neglect in Eq. ( 98 1 the contributions to the self- 



energy from intraband scattering current, e.g. via in- 
teraction with phonons, low energy photons (free car- 
rier absorption) or ionized impurities, since they cancel 
upon energy integration over the band. However, the 
current flow in a given band is not unaffected by the 



intraband scattering, since the carrier GFs are modi- 
fied by the action of the corresponding self-energies in 
the associated Dyson equations. 

Due to its generality, the above representation of 
the global current components in relation to the carrier 
self-energies has manifold applications, depending on 
the nature of the self-energies. First of all, it can be 
used to express the terminal current at contact a in 
terms of the contact self-energies, 



dE 



-tr 



E%<(E)G>(E)-EZ>(E)G<(E) 



(101) 



For an equilibrium contact characterized by a chemical 
potential /j, a and broadening function r R = i(Sa > — 
I7^ < ), the above current formula results in the Meir- 
Wingreen expression [9"2] 



dE 
2tt^ 



-tr 



r*(E){f, a (E)A(E)-G<(E)} 



(102) 



Another important application with relevance in the 
simulation of solar cells, especially for the so-called hot 
carrier devices, is the description of energy loss in elec- 
tronic transport via inelastic scattering processes such 
as coupling to phonons |29] , where the dissipated power 
is described in terms of the inelastic scattering current 
via En 



pdiss = _ 2 / dEJ ^ E ) E _ 



(103) 



But most importantly, the formalism provides the opto- 
electronic conservation laws in terms of the optical rates 
and associated currents via insertion of the electron- 



photon coupling self-energies in ( 98 ). On the other hand, 



the global radiative interband scattering rate equals the 
total optical rate R° pt , which, via the quantum statisti- 
cal average of the Poynting vector operator or directly 
from the explicit electron-photon self-energy terms in 
( 98 ) , can be expressed in terms of photon GFs and self- 
energies [53113], 

i?f = f d(MA / d 3 r f d 3 r'Y, p7V(r,r',M 

x ^<(r', r, M - (r, r', M^V, r , M 

(104) 

and which is related to the global electromagnetic power 
dissipation through 



p: 



■diss 



dui 
2^ 



hujR° pt (huj). 



(105) 
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The two terms in Eq. (104) provide microscopic expres- 



sions for absorbed and emitted photon flux, which may 
be regrouped to write the optical rate in terms of the 
net absorption, i.e. the absorption minus the stimulated 
emission, and the spontaneous emission, 



^L(M=^tr[^(M^ 7 ' < (^)], 
^ m (M=^tr[^<(/^)^(/^)], 



(106) 
(107) 



where the trace is over spatial coordinates and the carte- 
sian components of the vector potential, and we have 
introduced the spectral functions 



P 7 = £>7,> 



(108) 



for electron-hole pairs and for photons, respectively. 
The above expressions for -R° P * bs / em may be interpreted 
as global generation and recombination rates [HI]. The 
emission term provides a more general formulation of 
generalized Kirchhoff law |94j normally used in photo- 
voltaics for the photon flux emitted from an excited 
semiconductor. It allows for the inclusion of a realistic 
electronic dispersion via the polarization function #r, 
reflecting the effects of electron-electron, electron-hole 
and electron-phonon interactions and the non-equilibrium 
occupation of these states, as well as the consideration 
of the optical modes specific to the geometry of the 
system. 

In the electronic picture, the separation of the differ- 
ent contributions is not similarly straightforward, since 
the in- and outscattering self-energies always contain 
both an absorption and an emission component, e.g. 
interband absorption and intraband (free-carrier) emis- 
sion in the case of the electron inscattering self-energy. 
If the intraband terms are neglected, absorption and 
emission terms can be identified with the in- and outscat- 



tering contributions in ( 98 ) , where now the emission 



comprises both stimulated and spontaneous components. 

In the electronic expression for the optical rate, the 
GFs correspond to the final state of the scattering pro- 
cess, but the self-energy contains the GFs of the initial 
state only in the case of direct optical transitions. In 
the case of indirect, phonon assisted transitions such as 
those in indirect band gap semiconductors, the scatter- 
ing process proceeds via an intermediate, virtual state 
and the associated GFs and self-energies, as shown in 
Fig. [3] for the example of a simple three-band model of 
silicon [95] . 

Finally, the general formalism to treat scattering 
within the NEGF formalism can also be used to find 
an approach to the description of non-radiative recom- 
bination, e.g. via defects, which dominates in most pho- 
tovoltaic devices. While an exact microscopic treatment 
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Fig. 3 Self-consistent computation of Green's functions and 
scattering self-energies in silicon enabling the description of 
phonon-assisted indirect optical transitions |95 |. A similar 
procedure can be used for the simulation of non-radiative 
recombination via defects by replacing the virtual with a real 
midgap state. 



is possible only in special cases, a more phenomenolog- 
ical treatment proceeds along the lines of the general 
formalism introduced above. In order to describe non- 
radiative recombination via defects, the scattering in- 
duced renormalization of the GFs of carriers in localized 
as well as extended (bulk) states needs to be quantified 
by the corresponding self-energy expressions, which are 

k J 

S ^ k > E ^=J ^- h MZ /cap \KE)Gt{E T E), (109) 

where Gd is the GF of the defect or midgap state, Ai 
the matrix element of carrier capture or emission by 
the defect, and b = c,v for coupling of the defect to 
conduction and valence band states, respectively, and 
the upper sign corresponds to the former case. Inserting 
these self-energy terms into the general expression for 
the electron-capture rate yields 

i? c ^ = £ / ^.S> d (k,E)G<(k,E) 

k J 

k - 



J ^ h G>{E')S<{E')=R^ c 



(110) 



expressing the fact that the rate of scattering out of 
the extended conduction band states by carrier capture 
into defects equals the rate of scattering into the de- 
fect state. Using quasi-equilibrium approximation for 
the carrier GFs, the defect recombination rate derived 
by Shockley, Read and Hall [961197] is recovered. 

3 Numerical implementation 

For the numerical evaluation of the dynamical equa- 
tions Q, GFs and self-energies need to be represented 
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in a basis of finite dimension via a suitable expansion 
of the corresponding field operators, which converts the 
former in a linear system to be solved self-consistcntly 
with the diverse self-energy equations and the Poisson 



equation (75) for the electrostatic potential. Since the 



specific implementation depends strongly on the device 
configuration under consideration, some rather general 
issues arising in the numerical treatment of the NEGF 
formalism for solar cells will be discussed here, as well 
as the implications of commonly used approximations. 

3.1 Basis 

There are several aspects that may have an impact on 
the specific choice of basis for the carrier field oper- 
ators. First of all, there is the electronic structure of 
the material system: while direct gap semiconductors 
like the III-V materials commonly used in optoelec- 
tronic applications are reasonably well described by k p 
approaches, the situation is more demanding in indi- 
rect semiconductor such as silicon-based materials often 
used in photovoltaics, wherefore atomistic approaches 
such as empirical or ab-initio tight-binding are pre- 
ferred. For molecular devices like organic solar cells, 
ab-initio, quantum chemical methods such as density- 
functional theory are required for the electronic-structure 
calculation. There is an important difference between 
light-emitting optoelectronic devices and solar cells to 
consider at this point: while both transport and emis- 
sion take place in a narrow energetic range around the 
band edge, absorption of solar radiation is a broadband 
process. Hence, in principle, a full-zone description of 
the electronic structure is required to cover the photo- 
generation induced by the entire solar spectrum. 

Apart from the range of energy states described, 
the largest impact of the choice of basis concerns the 
spatial resolution and flexibility therein. While effective 
continuum approaches such as EMA or k • p lend them- 
selves to optimized adaptive spatial grids, allowing for 
coarse-grained treatment of spatially homogeneous re- 
gions such as bulk leads and high resolution close to 
heterointerface, the grid is fixed by the atomic lattice 
in the case of an atomistic basis, which precludes the 
description of extended systems, unless one resorts to 
multiscale approaches [55] . 

In some cases, a mode space basis may be preferred 
over a real space basis, i.e. in the case where only a 
few modes contribute to the device characteristics. This 
is for instance the case for the transverse modes in 
nanowire transport |18j . or in periodic superstructures 
[55] . However, in most cases, the mode spectrum in 
leads or contacts deviates strongly from the discrete 
mode spectrum of the low dimensional absorber, or, 



like e.g. in quantum well solar cells, there is a contin- 
uous transition from a discrete to a continuous mode 
spectrum in the energy range of interest. Indeed, the 
ability to treat on equal footing bound state to bound 
state and bound state to continuum transitions repre- 
sents a considerable advantage of real space approaches 



3.2 Energy discretization 

Among the numerical issues arising and approxima- 
tions imposed due to computational limitations, the 
most prominent ones for the photovoltaic applications 
in mind are inadequate energy mesh discretization and 
the non- locality of the interactions. The former problem 
arises whenever resonances appear in the transmission 
function of a structure, either due to localization effects 
leading to the formation of quasibound states or due to 
van Hove singularities in the density of states. The mesh 
should be refined around the resonances in order to 
avoid large errors in the energy integration and result- 
ing convergence problems. Different approaches have 
been developed to tackle the refinement issue, mainly 
adaptive algorithms [59,100] or methods for the reso- 
lution of the poles of the GFs giving rise to the reso- 
nances [S3] . In the case where a large number of poles 
are present, adaptive refinement according to the res- 
olution of the spectral output quantities is preferrable. 
However, in the presence of inelastic interactions such 
as coupling to broadband illumination, the situation 
may arise that any two energies are coupled, extend- 
ing the need for mesh refinement to the whole energy 
range and thus making the use of a coarse-grained start- 
ing mesh questionable. In any case of an irregular grid, 
care has to be taken to use interpolation schemes that 
guarantee current conservation. 



3.3 Solution of the linear system for non-local 
interactions 

The determination of the full GF provides information 
on propagation and correlations between any two sites 
of the dicvicc and thus includes all non-local phenom- 
ena. However, such a computation is very time and 
memory consuming, wherefore approximate schemes with 
a limited number of non-zero off-diagonal GF elements 
are used, which then allows for the application of fast 
computational routines such as the recursive GF algo- 
rithm [16] but at the same time reducing the range 
of non-locality taken into account by the simulation. 
While the consequences are sizable even for rather local- 
ized intraband scattering mechanisms such as coupling 
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to phonons |101|]102j , they are severe in the case of the 
interband scattering due to coupling to photons, since 
in the latter case, due to the comparably long wave 
length of the light, the non-locality is not restricted 
by the matrix element of the coupling itself, but by 
the density of states which couple [1001164] . Thus, non- 
locality is only limited in the case of optical transitions 
between localized states, as in weakly coupled quantum 
dot or quantum well absorbers. 

4 Applications 

In the following, the kind of information provided by 
the NEGF formalism shall be illustrated by applica- 
tion to the simulation of two different generic quantum 
photovoltaic devices, namely a standard III-V mate- 
rial based single quantum-well photodiode 66, 1031 1104) 
and a silicon based superlattice absorber |105j . While 
the former device embodies all relevant structural and 
mechanistic features of the multi-quantum well solar 
cell, the latter corresponds to a prototypical ingredient 
of third-generation multi-junction devices. The two de- 
vices are fundamentally different in several aspects of 
their working principles. First of all, in the QWSC, the 
confined states act only as absorbers, while transport 
takes place in the extended quasi-continuum states of 
the host material. In the superlattice device, both ab- 
sorption and transport proceed via the states of the 
coupled quantum wells. A second important difference 
is the role of scattering with phonons: while this interac- 
tion is a crucial aspect of the carrier escape mechanism 
in the QWSC, it greatly enhances transport in the su- 
perlattice at finite internal field due to restoration of 
resonant interwell charge transfer. It is thus instructive 
to have a closer look at the microscopic picture of charge 



carrier photogeneration and transport in these two dif- 
ferent devices, provided by the NEGF formalism in the 
form of spectral currents and rates. In both examples, 
coupling to equilibrium bulk phonons and isotropic free 
field photons is assumed. 

4.1 Carrier escape in direct gap single quantum well 
pin- diodes 

The escape of charge carriers from quantum wells under 
illumination and at finite terminal bias voltage is an ex- 
cellent example of a quantum photovoltaic process that 
can only be described properly by a quantum-kinetic 
approach, since neither the exact non-equilibrium car- 
rier distribution in the QW nor the phonon-assisted 
tunneling rate can be determined consistently in a semi- 
classical or ballistic non-local approximation. The model 
system used for the investigation of the carrier escape 
process, consisting of a p-i-n diode of high band gap 
material with a single quantum well formed by a slab 
of low band gap material allocated at the center of the 
intrinsic region, is shown in Fig. [4j The quantum well 
thickness is 18 monolayers (ML), with the ML width 
given by half the lattice constant of the semiconduc- 
tor material. The dopant concentration in the doped 
layers of 30 ML extension amounts to N^^ = 10 18 
cm~ 3 . The electronic structure is described by a sim- 
ple two-band tight-binding model [106] providing band 
gaps E\ = 1.77 eV and E™ = 1.42 eV and band offsets 
AE C = 0.2 eV, AE V = 0.15 eV. The parameters for 
electron-photon and electron-phonon interactions cor- 
respond to the GaAs-Alo.3Gao.7As material system (see 
Ref. jl03j for details) . Doping and band offsets result in 
a band diagram in growth direction as shown schemat- 
ically in Fig. [5] At the operating point of the solar cell, 
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Fig. 6 Spectral properties of charge carriers in the quantum well region of the device from Fig. [4] at contact Fermi-level 
separation of 1.3 V and under illumination with monochromatic light at 1.55 eV and intensity of 17.7 kW/m 2 : (a) Spatially 
resolved radiative generation and recombination spectrum, (b) spectral excess carrier density with phonon satellites, (c) electron 
in-scattering rate due to coupling to phonons, showing non-radiative charge transfer to higher states, (d) resulting current 
spectrum, with main contributions well above the energy of photogeneration |103 |. 



i.e. the device state where the product of current and 
voltage is maximum, the chemical potentials fj, n ,p gov- 
erning the injection of electrons and holes at the re- 
spective contacts are split by the corresponding bias 
voltage, but the occupation of the quantum well states 
is determined by the complex balance of the rates for 
photogeneration, radiative recombination, carrier cap- 
ture and carrier escape. Fig. [6] displays these rates to- 
gether with the resulting excess carrier concentration 
and current spectrum at a contact Fermi-level splitting 
of 1.3 V and under illumination with monochromatic 
light at 1.55 eV and high intensity of 17.7 kW/m 2 . The 
scattering rates are for electrons, a positive value thus 
corresponds to inscattering of electrons in the conduc- 
tion band and to outscattering of holes in the valence 
band. Due to the monochromatic illumination, the gen- 
eration spectrum is narrow, but the interaction with 
phonons leads to a clear spectral separation of absorp- 
tion and emission lines on the one hand (a), and to 
the formation of phonon satellites in the spectral den- 
sity of excess carriers on the other hand (b). Current 
flow is not centered on the energy of generation, but 
on much higher energy, as the wide triangular barri- 
ers prevent a direct tunneling escape from deep QW 
states (d). The energy resolved phonon scattering rate 
(c) clearly evidences the phonon mediated scattering of 
carriers out of the low energy states where they are gen- 
erated and into the less localized higher energy states 
with enhanced tunneling escape rate. The NEGF for- 
malism is thus able to provide a complete microscopic 
picture of the thermionic emission process at arbitrary 
operating conditions. 



4.2 Generation and transport in Si-SiO^ superlattice 
absorbers 

Owing to the large potential barriers for photoexcited 
charge carriers (>3 eV for both electrons and holes), 
transport in Si-SiO^ superlattice devices is restricted to 
quantum-confined superlattice states. Due to the finite 
number of wells and large built-in fields, the electronic 
spectrum can deviate considerably from the minibands 
of a regular superlattice. By including the coupling of 
electrons to both photons and phonons, the NEGF the- 
ory is able to provide a microscopic picture of indirect 
generation, carrier relaxation, and inter-well transport 
mechanisms beyond the ballistic regime. The model sys- 
tem used to investigate generation and transport in in- 
direct gap superlattice absorbers is shown schematically 
in Fig.[7^a). It consists of a set of four coupled quan- 
tum wells of 6 ML width with layers separated by oxide 
barriers of 3-ML thickness, embedded in the intrinsic 
region of a Si p-i-n diode. The thickness of both doped 
and spacer layers is 50 ML. The monolayer thickness is 
half the Si lattice constant. For such thin oxide layers, 
the confinement and thus the effective barrier height 
is reduced, and was chosen at 0.95 eV for both elec- 
tron and holes, which is still large enough to prevent 
thermionic emission. The carriers are described using 
a multiband effective mass model, accounting for the 
different conduction band valleys as well as light and 
heavy holes (details in |105j ). The doping density is 
Ni = 10 18 cm -3 for both electrons and holes. This com- 
position and doping leads to the band diagram shown in 
Fig.[7lb). In order to investigate the relaxation of photo- 
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Fig. 7 (a) Spatial structure and doping profile of the p-i(SL)-n model system. The doping level is Nd = 10 18 cm 3 for both 
electrons and holes, (b) band diagram of the p-i(SL)-n model system with the active quantum well absorber region | 105| . 



generated carriers, the superlattice region is illuminated 
with a photon energy of 1.65 eV, thus exceeding the ef- 
fective band gap, which at the chosen barrier thickness 
lies around 1.3 eV. At this photon energy, both first 
and second electronic minibands are populated, as can 
be seen in Fig. displaying the spatially and energeti- 
cally resolved photogeneration rate at an intensity of 
10 kW/m 2 and at vanishing bias voltage. Remarkably, 
as shown in Fig. [9j also the electronic current flows in 
both minibands, relaxation due to scattering is thus not 
fast enough to confine transport to the band edge. It is 
however obvious that transport of photocarriers across 
the superlattice is not ballistic, but strongly affected 
by the inelastic interactions and best described by the 
sequential tunneling picture. 



It remains to consider the indirect photogeneration 
process. To this end, the phonon-mediated intervalley 
in- and outscattering rates are displayed in Fig. [10] for 
both final electron state in the Si valleys close to the 
X-point and the virtual intermediate state at the r~ 
point. While in-and out-scattering are balanced in the 
dark (a) , there is a strong net charge transfer from r to 
X under illumination, and the inscattering rate reflects 
the spectral pattern of the photogeneration. 



5 Conclusions 

The high efficiency concepts of the next generation of 
photovoltaic devices pose new requirements on models 
for the theoretical description of their optoelectronic 
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Fig. 10 Spatially and energy resolved electron-phonon inter- 
valley scattering rate at short circuit conditions (a) in the 
dark and (b) under monochromatic illumination with energy 
E 1 = 1.65 eV and intensity I 1 = 10 kW/m 2 . 



characteristics. A theoretical framework suitable for this 
purpose is available in the NEGF formalism, if the quan- 
tum optics formulation is unified with the picture of 
dissipative quantum transport. The resulting theory is 
then able to provide insight into the modification of the 
central physical processes of photogeneration, charge 
separation and carrier extraction due to the use of semi- 
conductor nanostructures as absorber and conductor 
media. The general consideration all of the degrees of 
freedom on equal footing furthermore presents a suit- 
able starting point for consistent approximations in cases 
where some of the subsystems allow for a simplified 
treatment, as e.g. for bulk modes, equilibrium distribu- 
tions or negligible coupling parameters. While the gen- 
eral framework is set, many challenges remain for the 
application to realistic photovoltaic devices, both in ad- 
vanced but fundamental aspects such as excitonic con- 
tributions to the optical transitions and non-radiative 
recombination, as well as in the computational efficiency 
in treating extended system without excessive loss of 
accuracy. 
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